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Abstract 

The algebraic Bethe ansatz can be performed rather abstractly for whole 
classes of models sharing the same 7?-matrix, the only prerequisite being 
the existence of an appropriate pseudo vacuum state. Here we perform the 
algebraic Bethe ansatz for all models with 9x9, rational, gl(l |2)-invariant 
i?-matrix and all three possibilities of choosing the grading. Our Bethe 
ansatz solution applies, for instance, to the supersymmetric t-J model, the 
supersymmetric U model and a number of interesting impurity models. 
It may be extended to obtain the quantum transfer matrix spectrum for 
this class of models. The properties of a specific model enter the Bethe 
ansatz solution (i.e. the expression for the transfer matrix eigenvalue and 
the Bethe ansatz equations) through the three pseudo vacuum eigenvalues 
of the diagonal elements of the monodromy matrix which in this context 
are called the parameters of the model. 
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1 Introduction 

This work resumes previous work of one of the authors [1] where the algebraic Bethe 
ansatz for the gl(l |2) generalized model was constructed for the grading (+,—,—). In 
this article we address the two remaining cases (—,+,—) and (—,—,+) which turned 
out to be technically more involved, since the grading enters the auxiliary second level 
Bethe ansatz in a non-trivial way (see appendix iBl. 

Performing an algebraic Bethe ansatz calculation means to diagonalize the transfer 
matrix of a certain two-dimensional classical vertex model by purely algebraic means 
or (in a physicists language) by using only commutation relations between operators. 
If the transfer matrix has a Hamiltonian limit this is equivalent to diagonalizing the 
Hamiltonian along with its conserved currents. The Hamiltonian and the conserved 
currents are then usually generated by expanding the logarithm of the transfer matrix 
in the spectral parameter. 

The algebraic Bethe ansatz can be performed on a considerable level of abstraction 
and seems to depend only on the structure of the 7?-matrix of a given model and on the 
existence of a so-called pseudo vacuum or highest vector [2] on which the monodromy 
matrix acts as an (upper) triangular matrix. This idea was first of all worked out for 
the models falling into the same class as the XXZ spin chain [3] and turned out to 
be useful in calculating the norm [4] and certain matrix elements [5] of Bethe ansatz 
states. In [3,4] V. E. Korepin introduced the notion of a 'generalized model' whose 
'representation' is given by the action of the diagonal elements of the monodromy 
matrix on the pseudo vacuum. He assumed the vacuum eigenvalues, say ai{X) and 
a2{'k), to be arbitrary and called them (functional) parameters of the model. Later [6,7] 
V. O. Tarasov refined and basically confirmed Korepin's concept. 

The algebraic Bethe ansatz for the generalized model associated with the 7?-matrix 
of the XXZ model is of the same structure as for the fundamental XXZ spin chain. 
Such simple relation holds no longer for models of 'higher rank' which require a nested 
Bethe ansatz. The simplest models which allow for a nested algebraic Bethe ansatz are 
the models with gl(n) invariant 7?-matrix [8,9]. Considering the fundamental represen- 
tations of these models one observes that not only the monodromy matrix elements be- 
low the diagonal annihilate the pseudo vacuum, but additional zeros appear above the 
diagonal [10]. This fact simplifies the algebraic Bethe ansatz for the fundamental rep- 
resentation as compared to the more general case, where the action of all the elements 
of the monodromy matrix above the diagonal is non-trivial. For the solution of this 
more general case a new concept, the vacuum subspace, was introduced by Kulish and 
Reshetikhin [11]. This new concept made it possible to perform the algebraic Bethe 
ansatz for the models with gl(n) invariant 7?-matrix on the same level of generality as 
in the gl(2) case (corresponding to the 7?-matrix of the XXX spin chain of spin ^). The 
resulting eigenvalue of the transfer matrix and the Bethe ansatz equations depend on 
n functional parameters ai(?i), . . . ,a„(?i), which, together with the triangular action of 
the monodromy matrix on the pseudo vacuum, define the gl(n) generalized model [12]. 
Considering the parameters ai{X),a2{'k),a^{X) as free Reshetikhin derived the norm 
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formula for the gl(3) case [12]. 

To our knowledge there was no more activity in the direction of constructing alge- 
braic Bethe ansatz solutions of generalized models beyond the above mentioned work 
of Kulish and Reshetikhin. This may partially be related to the general difficulties in 
generalizing the algebraic Bethe ansatz beyond gl(n) (see e.g. [13]). For the models 
with gl(l|2)-invariant i?-matrix (e.g. [14-22]^), which (together with its anisotropic 
generalizations) are of particular interest for application in solid state physics, many 
algebraic Bethe ansatz solutions were constructed [20,21,23-29] which can all (except 
for [26]) be interpreted as certain realizations of the solution obtained in [1] and are 
extended here to the two gradings (—,+,—) and (—,—,+) not treated in [1]. 

Future applications of our work may be the quantum transfer matrix approach to the 
thermodynamics of the models with gl(l|2) invariant i?-matrix [30,31], the calculation 
of norms and matrix elements and, possibly (see [1]), the algebraic Bethe ansatz for 
the Hubbard model. 

Our article is organized as follows. Section 2 contains the basic definitions relating 
to the gl(l |2) generalized model with some material shifted to appendix|X| In section 3 
we discuss how the grading and the Yang-Baxter algebra change under permutations 
of the basis vectors in auxiliary space. We shall see how the grading may change with 
the change of the highest vector. Section 4 contains our main result, the formulae 
for the algebraic Bethe ansatz solution of the generalized model for the three possible 
gradings. For the grading (+,—,—) a derivation was presented in [1]. The proof for the 
remaining two cases is sketched in appendix |Bj Section 5 contains examples of how 
to apply the formulae of section 4. We mainly reconsider the well known example of 
the fundamental supersymmetric t-J model [24,25] and clarify the connection between 
the possible choices of pseudo vacua and possible gradings. We also obtained a simple 
proof of the equivalence of the different Bethe ansatz solutions of the supersymmetric 
t-J model (which exceeds the one in [24], since we also show that the eigenvalues 
are identical). This proof [32] will be published separately. Section 6 contains our 
conclusion and a discussion of future perspectives. 

2 The gl(l 1 2) generalized model 

We begin by specifying the class of models we are going to consider. As was explained 
in the introduction this class is determined by its i?-matrix and by the existence of a 
pseudo vacuum. Here the 7?-matrix has matrix elements of the form 

i?g(X) =a(X)(-l)^(«)^W5^5j + M^)5^5j, (1) 

a, (3, Y, 5 = 1,2, 3. This i?-matrix is based on graded permutations [15]. It is contained 
in the early list of Kulish and Sklyanin [8]. The algebraic Bethe ansatz for the funda- 
mental model was constructed by Kulish in 1985 [23] along with the general gl(m|n) 
case. 

a more thorough discussion see the introduction of [1]. 
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The i?-matrix (HI) is characterized by two rational, complex valued functions 

a{X) = ^, b{X) = ^ (2) 
A + ic A + ic 

depending on a spectral parameter X E C and a coupling c EC It further depends on 
the grading, /) : {1,2,3}^Z2. We shall consider three cases:* 

(i) pil)=0, p{2)=pi3) = l, 

(ii) p(2)=0, p(3)=p(l) = l, (3) 

(iii) ;?(3)=0, p{l)=p{2) = l. 

In order to refer to these different cases we introduce a vector index g which is g = 
{gi,g2,g3) = in the first case, g = (—,+,—) in the second case, and g = 

(—,—,+) in the third case. We shall say 'the grading is g\ and we shall write Rg 
instead of R. 

The matrix RgCk) solves the Yang-Baxter equation and obviously satisfies the com- 
patibility condition [8] 

= (-1)M«)+/^(P)+/^W+MS)^^«J(X) . (4) 

In order to introduce the notion of the graded Yang-Baxter algebra we shall further 
need the matrix RgiX) defined by switching the row indices of RgCk), 

^s'^li'^) = ^g^^C^) ■ (5) 

The graded Yang-Baxter algebra with 7?-matrix R{X) is the graded, associative al- 
gebra (with unity) generated by the elements T^{X), a, (3 = 1,2,3, of the so-called 
monodromy matrix modulo the relations 

R{l-lj){T{l)(^gl{p)) = {'T{^i)^gT{l))R{l-ij). (6) 

We shall assume that the elements of the monodromy matrix are of definite parity, 
7r(Tp°'(?i)) = p{a) + p(|3). The symbol ®g denotes the super tensor product associated 
with the grading g. For a definition see appendix A. 

We now define the set of all models related to the i?-matrix ([T]) and solvable by 
algebraic Bethe ansatz as 'the gl(l|2) generalized model': By definition the gl(l|2) 
generalized model is the set of all (linear) representations of the graded Yang-Baxter 
algebra © having a highest vector (or pseudo vacuum) Q.. The highest vector Q is a 
vector on which the monodromy matrix T (X) acts as an upper triangular matrix and 
which is an eigenvector of its diagonal elements: 

(V) 

T^{X)^ = Q, for a > (3. 



*We comment on the remaining possibilities of choosing the grading below in section|3] 
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The eigenvalues cijiX), j = 1,2,3 of the diagonal elements of 'T{X}, are called the 
parameters of the generalized model. These parameters characterize the representation 
in a similar manner as the highest weight in a highest weight representation of a Lie 
algebra. 

Let us denote the representation space of a given representation of the generalized 
model hy ^ . It is clear from the quadratic commutation relations contained in the 
graded Yang-Baxter algebra © and from ^ that we may assume ^ to be spanned by 
all vectors of the form 

4.(^1 , . . . , ^;v) = r^^; ill)... T^l' . (8) 

where ak < ^k, k = I, . . . ,N. This assumption is at least sensible for a finite dimen- 
sional representation space i?/ . 

The super trace of the monodromy matrix 

t{X) = {-ir^^^'r^{X) = str,{'r{X)). (9) 

is called the transfer matrix of the generalized model. Since R{X) is invertible for 
generic values of X G C, we conclude from © and Q that the transfer matrix satisfies 

[t{l),t{p)] = (10) 

for all generic X,^ EC. It follows that t(k) and t{ij) have a common system of eigen- 
functions, which means that the eigenvectors of t{X) are independent of the spectral 
parameter X. 

The task to be performed below of the algebraic Bethe ansatz for the generalized 
model is to diagonalize t{X), i.e., to solve the eigenvalue problem 

t{X)^ = A{X)^. (11) 

It is a remarkable fact that this task can be accomplished by solely using the graded 
Yang-Baxter algebra ® and the properties Q of the highest vector Q.. In particular, 
it is not necessary to require that T^{X)Q. = as in case of the fundamental graded 
representation, which corresponds to the supersymmetric t-J model. 

3 Variation of the grading 

Before presenting our results for the algebraic Bethe ansatz we would like to explain 
why we may restrict ourselves to upper triangular action in our definition Q of the 
highest vector Q. and why we consider only the three gradings shown in equation Q. 
For the former purpose we first of all introduce the natural action of the symmetric 
group 6^ on row vectors x = (xi,a:2,X3), setting 



(12) 
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Figure 1: Change of the monodromy matrix action on the highest vector Q. under 
permutations of the basis vectors in auxiliary space. 

for all o E This defines a faithful representation of &^ which is orthogonal with 
respect to the usual Euclidian scalar product (x,y) = xiyi +X2y2 +^3^3- 

Denoting the transposed of a by & we obtain the transformation properties of the 
i?-matrix under permutations directly from its definition ([T]), 

{o'®o')Rg{l){o0(5)=Rgail). (13) 

Similarly, the graded Yang-Baxter algebra ^ is easily seen to transform as 

Rga{X-iu){a'T{X)o®g^o'T{ij)a) = {a'T{ij)o(»g^o''r{X)o)Rg^{X-fj), (14) 

while it follows from Q that 

str^(T(;i)) = strg,,(a'T(;i)a), (15) 

which expresses the invariance of the transfer matrix with respect to permutations. 

In figure[T]we show the action of the transformed monodromy matrix o'T {TCjo on 
a highest vector Q., when o runs successively through all permutations in generated 
by the transpositions of nearest neighbours and 1123. We see that if, for a given 
grading g, a monodromy matrix T (k) realizes one of the six patterns in figure [l] by 
acting on some vector Q. then there is a permutation o E &^ such that &T {X)(5 acts as 
an upper triangular matrix on Q.. The corresponding grading changes from g to go. 

Let us consider an example. Take g = (— , — , +) and T (X) and Q. such that 

/fli 0\ 

T{'k)a= * a2 * (16) 
V* asj 

This is the pattern in the middle of the second row in figure [T] Thus, 

(a2 * *\ 
as * \n, (17) 
ai/ 
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and the grading changes to (-, -, +)ni2n23 = (-, +, -)• 

Combining the six patterns in figure [T] with the three gradings g = (+,—,—), 
(—,+,—), (—,—,+) we obtain 18 cases which by application of permutations all re- 
duce back to three, e.g. upper triangular action with three different gradings. The Bethe 
ansatz solutions of the transfer matrix eigenvalue problem (fTTT) for these three cases 
will be presented in the next section. 

Note that it may happen that there are several vectors Q.i,Q.2,... which for given 
monodromy matrix and grading generate several of the patterns in figure ([T]). Then 
there are several equivalent but differently looking Bethe ansatz solutions of the trans- 
fer matrix eigenvalue problem. This is, for instance, the case for the supersymmetric 
t-J model as was observed in [24]. We will come back to this phenomenon in our 
example section (see section|5l). 

How about the other possible gradings? There are eight cases altogether. Three 
cases are listed in ©. The case g = (+, +, +) was treated by Kulish and Reshetikhin 
[11]. The four remaining cases are related to our work or the work of Kulish and 
Reshetikhin by a switch of sign of g (e.g. — , — ) (— , +, +)). It is easy to see 
that this switch modifies the Yang-Baxter algebra Q only trivially: We introduce a 
diagonal matrix G = diag(^i,^2,^3) and the 3x3 unit matrix 73. Then a similarity 
transformation with G0 13 transforms the graded Yang-Baxter algebra Q with R- 
matrix Rg into 

R-gil-lj){'T{-l)®^gr{-iu)) = (18) 

Note that the expression for the parity of the monodromy matrix elements, n{T^{X)) — 
p{a) + /'(P), is invariant under a change of the sign of g, since it corresponds to re- 
placing p{a) by p{a) + 1. Thus, every representation with parameters aj{X} of the 
graded Yang-Baxter algebra Q with grading g and i?-matrix RgQCj is equivalent to a 
representation with parameters aj{—X) of Q with grading —g and i?-matrix R-g{X). 
Consequently the Bethe ansatz solutions of the generalized model for the remaining 
gradings are obtained from the solutions in the following section by switching the sign 
in the argument of a(X) and the overall signs of the transfer matrix eigenvalues (or by 
performing similar manipulations in the corresponding equations in [1 1]). 

4 The algebraic Bethe ansatz solution 

In this section we list the transfer matrix eigenvalues and the corresponding eigen- 
vectors. For g = (-(-,—,—) they were obtained in [1]. The derivations for the two 
remaining cases g = (— , — ) and g = (— , — , +) are presented in appendix iBl 

Recall that the functions ai, a2 and are the (functional) parameters of the model 
and that a{X) = X/ (?i + ic). The different transfer matrix eigenvalues Ag(A.) (equations 
(fT^ . (ETT) . (1^ below) are distinguished for a given grading by specifying two sets of 
Bethe roots {kj}^^^^ and {^k]f=\ which have to be calculated from two coupled sets 
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of Bethe ansatz equations (see (l20b . (l22t . (l24l) below). 



A(+__)(X) =ai(X)n 



W 1 M 1 Ml 



^-nV-^' J = 1,.-,A^, (20a) 

l^k 



N -I M Ml 

y a (X - X;) [nk - A) a [nk - A) 



^ = n^^' i = l,.-,A^, (22a) 
^ = f[^-^rT. *=1.-.M. (22b) 



A(__+)(X)--ai(X)n- ^ 



iV 1 M 1 Ml 

• «2 (^) n ^0-^ n — ^ + «3 (^) n -n — v ' ^^s) 



aiiXj) fJya{Xi-Xj)lJ^a{Xj-iJk)' 

¥j 

^ = n^A^^ ^=l,...,M. (24b) 

These three sets of expressions for the eigenvalues and Bethe ansatz equations depend 
on the grading g = (^1,^2,^3) in a characteristic way which allows us two write them 
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all in one (for a similarly compact expression for the (^-deformed) fundamental model 
see [33]): 



N 



A' J M J M ^ 

1=1 «U2(A-Xj)) ^^-^ a{g2{pk-/^)) tJi ^iSii^-m) 



a,iXj)_^ aig,ih-lj)) f. 1 y=i,...,A^, (26a) 

_y^ a{g3iiJk-IJi)) A 1 j^^-^ ^ ^26b) 

i^k 

Our notation means, for instance, that by specifying gi = —I, g2 = I, g3 = —I equa- 
tions (l25b . (l26b turn into (l2Tt . (l22t corresponding to the grading g = (—,+,—). 

Describing the corresponding eigenvectors requires more effort, since we will have 
to introduce several notions related to the 'second Bethe ansatz' in the nested Bethe 
ansatz calculation that led to the above expressions for the eigenvalues. The eigenvec- 
tors are obtained by acting with certain linear combinations of products of monodromy 
matrix elements on the highest vector Q.. They are of the form 



4>^(Xi,...X^;^i,...,A/m)= £ bI{1,)...BI{InW^---'^{iuu...,ium)^. (27) 

'lvi'A'=l 

Here the 5y , j = 1,2, may be thought of as components of row vectors 

B(+-\l) = {r,\l),T,\l)), (28a) 
5(-+-)(X) = (T3i(X),T2i(^)), (28b) 
5(-+)(X) = (T2^(X),'r3i(X)). (28c) 

(Pg is defined in terms of an auxiliary monodromy matrix which is a product of two 
2 X 2-matrices 
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The factor 'Dg{X) basically contains elements of the monodromy matrix T {X), 



^^--)(^^^l.3^j£)- '^^SJ ''''' 
( T32(X)(0^)®n 

The factor T^'(?i) is the monodromy matrix of an auxiliary 'spin problem', 

T.KX) = 4(kM-'k)...L((k,-X), (31) 

carrying an induced grading g' , which is (+, +) for g = (+, — , — ) and (+, — ) for the 
remaining two cases g = (—,+,—) and g = (—,—,+). The corresponding elementary 
L-matrices are^ 

L^^^\l) = a{X)h + b{X) {'^\ '^^^ . (32) 

and 




Lf'\V)=a[V)h + h[V)l'\ _'\\. (33) 

The ej^ are the canonical basis elements of (End(C2)) ®^ (viewed as a graded algebra) 

introduced in appendix|Xl They depend on the induced grading g'. The L-matrix -^^j^^^ 

is the L-matrix of the XXX spin chain of spin 5, whereas Lj^ ^ generates the gl(l 1 1) 
invariant model of free spinless fermions. 

Inserting (l30b - (l33t into (l29t we have defined the auxiliary monodromy matrix 
T (X) for the three different gradings under consideration. The 2^-column vector cp^ 
in equation (l27t is constructed by acting with matrix elements of T (A,) on appropriate 
auxiliary states. For g = (+, — , — ) and g = (— , — , +) we define 

(Pg(^i,...,A/M) =5(a/i)...5(a/m)(J)®^, (34) 
for g = (—,+,—) an appropriate definition is 

(P^(^i,...,A/m) =C{pi)...C{pm)(D'^^ ■ (35) 

A derivation of the eigenvectors and the corresponding eigenvalues for g = (—,+, —) 
and g = (—,—,+) is presented in appendix |Bl The proof for g = (+,—,—) can be 
found in [1]. 

^Note that, as compared to [1], we have changed the definition of the L-matrix Lj^^'^ into the 



equivalent form yj] 
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5 Examples 

Numerous examples of systems for which the Bethe ansatz solution of the previous 
section applies can be constructed starting from the observation [23] that 

Lj^il)=a{l^ + b{l){~ir^''^P^^^Ej^ (36) 

with a{X) , bCk) defined in Q is a representation of the graded Yang-Baxter algebra Q 
if Ej^ is a representation of gl(l |2) of parity compatible with the grading, 

[^;>^4]± = ^A^m - (-i)(M«)+MP))(MY)+.(8))5y£^.a) _ (37) 

Here [■, -ji denotes the superbracket (see appendix IXb. 

For the construction of models of fermions on one-dimensional lattices (which is 
our personal concern with the Bethe ansatz presented in the previous chapter) one may 
generally utilize (l36b in a way that involves three logically separate steps: 

(i) Take a representation of gl(l |2) in End(C") (eq. ^ with ; = k, see [34, 35]). 

(ii) Embed it into (End(C"))®^ in such a way that the grading ( (l37t for j ^ k) is 
respected [36]. 

(iii) Introduce Fermi operators [36,37]. 

Since we merely want to illustrate our Bethe ansatz solution of the previous sections 
we shall take the three steps in one in the examples considered below. 

The most elementary example is, of course, the supersymmetric t-J model [16,24] 
which is the fundamental model associated with the i?-matrix dl]). The supersymmetric 
t-J model is a model of electrons on a lattice. In order to be able to write down the L- 
matrix and the Hamiltonian in a familiar way we introduce canonically anticommuting 
creation and annihilation operators cj ^, Ck,b where the indices j,k = 1, . . . ,L refer to 
the lattice sites, and a, =t, i are spin indices. 

Due to the canonical anticommutation relations the elements (Xj)^, a, [3 = 1,2,3, 
of the matrix 



:i-«7i)(l-«7T) (l-"ji)CjT C7i(l-"iT)\ 



(38) 



with n j I = Cj |CjVf , Hj i^ = c ■ ^c^ ^ form a complete set of 'projection operators' on 
the space of states locally spanned by the basis vectors |0), cj||0), cj||0). Double 
occupancy of lattice sites is forbidden on this space. Setting X^P = (X;)p, a, |3 = 1, 2, 3, 



we find 



XjlXk^ = (_l)(Kcx)+p(|3))(p(Y)+P(8))x,5x,P , for j ^ k, (39b) 
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where p{l) = 0, p{2) = p{3) = 1. It follows that the operators Xj^ satisfy equation 
(l37t . The linear combination Xj'^ = 1 —nj-^nj^ projects the local space of lattice elec- 
trons onto the space from which double occupancy is excluded. The corresponding 
global projection operator is 



It will be needed below. 

We conclude with dSTJ) and ^ that the L-matrix 



(40) 



Lj{'k)=a{'k)l3+b{X) 



^J2 



-Y-2 



(41) 



is a representation of the graded Yang-Baxter algebra Q with grading (-f, — , — ). This 
representation has been termed fundamental graded representation in [36]. The Fock 
vacuum |0) satisfying Cy,a|0) =Ofor j= 1, . . . ,L; a,Z7 =t, | is clearly a possible highest 
vector for Lj{X), 



Lj{l)\0) 



b{X)Xj\ b{X)Xj^^ 
a{X) 
a{X) I 



|0). 



(42) 



It turns out that the matrix LjQC) generates the supersymmetric t-J model at a 
single site. The corresponding monodromy matrix of the L-site model is 



Its action on the Fock vacuum follows from (l42b as 

/l Bx{V) B2{1)\ 
a^{X) 



(43) 



Til)\0) 







|0). 



(44) 



Thus, T (k) is a representation of the graded Yang-Baxter algebra ®, and |0) is a 
highest vector satisfying ©. It follows that our general formulae (fT9b . (l20b apply with 
functional parameters which can be read off from (l44b : 

ai{X) = l, a2{X)=a3{X)=a^{X). (45) 

This way we have recovered equations (3.47), (3.48) and (3.50) of [24]. 
Note that the Hamiltonian of the supersymmetric t-J model is 



// = -icaxln{(str('r(0))) 'str('r(X))} 



A,=0 



(46) 
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Because it acts on the restricted space of electronic states, where no lattice site is 
doubly occupied, we may replace it with (see [36, 37]) 



where we inserted the usual definitions Sj = j 'Ea.h=tl a^ah'^j b local spin opera- 
tors in terms of Pauli matrices o*^ and introduced the local particle number operators 

Clearly the Fock vacuum is also a highest vector for 11231' {X)Yl23, since under a 
permutation of the second and third basis vector in auxiliary space Lj{X) |0) transforms 
into 

/l b{X)Xj\ b{X)Xj}\ 
aiX) 
^0 a{X) I 



n23i:,(^)n23|0) 



|0). (48) 



which is of upper triangular form. The monodromy matrices T (?i) and 11231' (X)n23 
carry the same grading and, as for arbitrary permutations, lead to the same transfer 
matrix (see equation dTSt ). Since the grading as well as the parameters ai, ^2 and a-^ 
are identical for T (?i) and 11231' (?i)n23, both monodromy matrices lead to the same 
form ([T9b of the transfer matrix eigenvalue and to the same Bethe ansatz equations 
(l20b . Note, however, that the Bethe ansatz eigenvectors (l27t are different, because 
B^^ ) (?i) is changed to B^^ ^ which is equivalent to a spin flip transformation. 
Indeed, applying the spin flip transformation 

7W = n(l-(ct^-c]j(c,,,-c,J) (49) 

to the elementary L-matrix we obtain 

7(-')i:,(^)(/(^))^ = n23i:,(^)n23 (50) 

which implies the invariance of the transfer matrix with respect to spin flips. 

It is a well known fact that there are three alternative sets of Bethe ansatz equations 
for the supersymmetric t-J model [24,25]. Let us see how this comes out in our general 
formalism. We shall consider the monodromy matrix action on the two states 

\\)=X^2...Xx\\^) = cl^y..c\^\^), (51a) 
\\)=XjX...X,\\<d) = cl,...c\M. (51b) 



14 



OG 6^ 




§ 




a2{X} 


am 


id 


|0) 







_) 


1 


a\X) 




ni2 


IT) 




+, 






a\X) 


a^{X) 


ni2n23 


li) 




+, 




aH-X) 


a\X) 


a^{X) 


111211231112 


U) 






+) 


a^{X) 
aH-X) 


a\X) 




n23ni2 


IT) 






+) 


a^(X) 
aH-X) 


a\X) 


a^{X) 


n23 


|0) 








1 




a^{X) 



Figure 2: Different Bethe ansatz solutions for supersymmetric t-J model. 



Calculating first of all the action of the L-matrix on these states we obtain 



/ ci{X) 



i^;(^)IT) 



^;Wli) 







\ 



b{X)Xj\ a{X)-b{X) -b{X)Xjl 

\ a{X) j 

I a{X) ^ \ 

a{X) 

\b{X)Xj] -b{X)Xjl a{X)~b{X)j 



T). 



i)- 



(52) 



(53) 



Comparing the patterns of zeros on the right hand side of these equations to figure[T]we 
see that | t ) is a highest vector for ni2T (?i)ni2 and for 112311121 (?i)ni2n23, while 
I I) is a highest vector for ni2n23T (?i)n23ni2 and for ni2n23ni2'T (?i)ni2n23ni2. 
Together with the possibilities already covered by using |0) as a highest vector we ob- 
tain all six cases of figure [l] albeit with different grading. The situation is summarized 
in tabular 121 

Taking the data from the tabular and inserting it into (l25t - (l27t we obtain the dif- 
ferent possible Bethe ansatz solutions of the supersymmetric t-J model. The expres- 
sions for the eigenvalues and Bethe ansatz equations are in agreement with the results 
of [24]. Because of space limitations we do not repeat those results here. In [24] 
and also in [25] the authors avoided writing explicit expressions for the L-matrix in 
terms of Fermi operators. Therefore they could not see the correspondence between 
the possible choices of pseudo vacua |0), | T )> I i ) the possible gradings. As we 
can learn from tabular El this correspondence is not unique. For the highest vector | ] ) 
the Bethe ansatz can be realized with grading (—,+,—) or (—,—,+), respectively. A 
similar statement holds for | |) . By way of contrast, the Bethe ansatz equations and 
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the expression for the transfer matrix eigenvalue are uniquely fixed if we choose |0) 
as the highest vector. Still, as was observed above, the eigenvectors can be realized in 
two different ways. 

The algebraic Bethe ansatz for the supersymmetric t-J model is rather special as 
compared to the general case. This is due to the fact that ajQO) = a3{X) (see tabular|2l) 
and that an additional zero appears above the diagonal in the action of the monodromy 
matrix on the highest vector. Consider, for instance, the case Q.= |0), o = id. Then 
the monodromy matrix action on Q. is given by (EH). Making contact with the notation 

of (QUI) we see that )(^)|0) = a^{X)l2\0). Because of this trivial action we may 

drop the factor 'Dg{X) on the right hand side of (12^ and the eigenvectors dTTl are 
written only in terms of Tr^ (X) , T^^ (X) and the auxiliary monodromy matrix ■ 
(see (EB, (El)). 

More examples are obtained by inserting other representations of (l37t into (l36t . 
For the grading g = (+,—,—) we may, for instance, take the four-dimensional repre- 
sentation 



^ ch^(a)— Hj— ch(a)c| — e "cjnj ch(a)cj — e "cJnA 
ch(a) — e^"c|n| cjc-]- 
\ch(a)q-e-«C|n| cjq J 



(54) 



which depends on a free parameter a G M. Note however, that it requires more effort 
[2] to associate a physically relevant model having a local Hamiltonian with higher 
dimensional representations. This is a subject that exceeds the scope of this work. 



6 Conclusion 

We have completed the task, begun in [1], of constructing the algebraic Bethe ansatz 
for the gl(l|2) generalized model. In this work the gradings g = (—,+,—) and g = 
(— , — , +) have been treated. As we hope to have convinced the reader in sections 3 
and 5, a complete understanding requires to consider the three gradings together. 

We hope our work will prove to be useful in future constructions of exact solutions 
of models with 7?-matrix ([T]). In first place we think of novel impurity models and 
of possible applications to Yang's model of electrons interacting via delta function 
potential and, maybe, to the Hubbard model (see discussion in [1]). Other applications 
may be the calculation of quantum transfer matrix eigenvalues (the quantum transfer 
matrix has a staggered pseudo vacuum on which the monodromy matrix acts 'without 
producing additional zeros above the diagonal') and the calculation of norms of Bethe 
ansatz eigenstates (see [1]). 
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Appendix A: Graded algebras 

In this appendix we shall recall the basic concepts of graded vector spaces and 
graded associative algebras. In the context of the quantum inverse scattering method 
these concepts were first utilized by Kulish and Sklyanin [8, 23]. 

Graded vector spaces are vector spaces equipped with a notion of odd and even, 
that allows us to treat fermions within the formalism of the quantum inverse scattering 
method (see [36, 37]). Let us consider a finite dimensional vector space V, which is 
the direct sum of two subspaces, V — Vo®Vi, dimVb = w, dimVi = n. We shall call 
vo G Vq even and vi e Vi odd. The subspaces Vq and V\ are called the homogeneous 
components of V. The parity 71 is a function Vi — > Z2 defined on the homogeneous 
components of V, 

%{vi) = i, / = 0,1, VieVi. (A.1) 

The vector space V endowed with this structure is called a graded vector space or super 
space. 

Let J? be an associative algebra (with unity), which is graded as a vector space. 
Suppose X, y e are homogeneous. If the product XY is homogeneous with parity 

Tl{XY) = %{X) + 7t(F) , (A.2) 

then R is called a graded associative algebra [8]. 

For any two homogeneous elements X, 7 e J? let us define the super-bracket 

[Z,F]±=XF-(-lfTO%X, (A.3) 

and let us extend this definition linearly in both of its arguments to all elements of si . 

Let p : { 1 , . . . , n} — >^ Z2. The set of all n x n matrices A, 5, C, . . . with entries in Jl , 
such that 7i(Ap) = 7t(Bp) = 7i(Cp ) = • • • = p(a) + p{^) is an associative algebra, say 

Mat(j?,n), since %{A^By) = p{a) + p{y). For A, 5 e Mat(jl,n) we define the graded 
tensor product (or super tensor product) 

(A®^5)^J= (-l)(^(«)+^(P))^(YU«5j. (A.4) 

The graded tensor product is associative. For matrices A,fi,C,D G Mat(jl,n) with 
[5p,C^]± = Owe find 

(A 0g B) (C (8)g D) = AC (8)g ED . (A.5) 

Our chief example of a graded associative algebra is the algebra Jl — (End(y))®^, 
where V is a graded vector space as introduced above. acquires the structure 
of a graded algebra in the following way: We fix a basis {ei, . . . of defi- 

nite parity and define p : {1, . . .m + n} ^ Z2 by setting p{a) = %{ea). Then the set 
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{4eEnd(y)|a,|3= l,...m + n} with 

^a^Y = ^y^a is a basis of End(y), and the tensor 
products ea\®---'^ span the vector space R = (End(y))®^. Hence, the definition 

7i(4'i ® ■ ■ ■ ® efc) = p(ai) + + ■ ■ ■ + piat) + p{^l) (A.6) 

induces a grading on Jl regarded as a vector space. It is easy to see that an element 
X = ■' j^^^aj ® ■ • ■ ® G ^ is homogeneous with parity 7l(X) if and only if 

^_i)LUp(^j)+p(^j))x^;--g^ = i-\)''^''^x^;::g' . (aj) 

But the latter equation implies that two homogenous elements X and Y satisfy equation 
(IA.2I) . and is a graded algebra. 

The following definition of 'graded local projection operators' [36] provides a stan- 
dard basis of the graded associative algebra which is of crucial importance in con- 
structing solvable lattice models, 

,^.P = (_l)(p(a)+p(P))l£;p(Y.),Yi ^,P^ (A.8) 

Here is the {m + n) x (m + n) unit matrix, and summation over double tensor 
indices (i.e., over yi, . . . .Jj-i) is implied. The index j on the left hand side of (IA.8I ) 
is called the site index. A simple consequence of the definition (IA.8I ) for 7 7^ ^ are the 
commutation relations 

^jWj = (-l)(M«)+/'(P))(/'«+M5))^^5^^.p _ 
It further follows from equation (IA.8I) that eyP is homogeneous with parity 

%{ejl)=p{a)+p{^). (A.IO) 



Hence, equation (IA.9I ) says that odd matrices with different site indices mutually an- 
ticommute, whereas even matrices commute with each other as well as with the odd 
matrices. For products of matrices ej^ which are acting on the same site (IA.8I) implies 
the projection property 

^7tiy = Sy^a- (A.ll) 
Using the super-bracket (IA.3I) . equations (IA.9I ) and (lA.l II) can be combined into 

[ele,'yh = (5?./„ - (-1) W«)+^(P))(^«+^(S))6S.,P) . (A.12) 

The right hand side of the latter equation with j = k gives the structure constants of the 
Lie super algebra g\(m\n) with respect to the basis {sj^}- 

Since any m + n-dimensional vector space over the complex numbers is isomorphic 
to C'"+", we may simply set V = C™+". We may further assume that our homogeneous 
basis {ea E C^+'^ja = 1, . . . ,m + n} is canonical, i.e., we may represent the vector ea 
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by a column vector having the only non-zero entry +1 in row a. Our basic matrices 
are then (m + n) x (m + n) -matrices with a single non-zero entry -i-l in row a and 
column p. 

The definition (IA.8I) generalizes the notion of the Jordan- Wigner transformation to 
systems with higher spin (see [37]). As with the Jordan-Wigner transformation another 
consistent definition of the graded local projection operators, also leading to (IA.9I ) and 
(lA.llI) . is obtained by placing the factors ( — l)(iH«)+p(P))p(Y<;)g^^ the tensor product 

on the right hand side of (IA.8I) behind rather than in front of e^- This alternative 
convention was used in [36, 37]. Here we use (IA.8I ) instead as it naturally appears 
in the derivation of the algebraic Bethe ansatz for the gl(l |2) generalized model with 
gradings g = {—,+, — ) and g = (— , — , +) (see appendix|Bl). 

Appendix B: Derivation of the eigenvectors and eigenvalues of the 
transfer matrix for g = (—,+,—) and g = (—,—,+) 

It is most convenient to start with the case g = (—,—,+), for which the Yang- 
Baxter algebra has a simple block structure, and to obtain the case g = (— , -f, — ) af- 
terwards. In fact, it is equivalent to the case ^ = (—,—, + ) but with a transfer matrix 
acting on the pseudo vacuum Q. as 

(ai{X} * * \ 

asil) a. (B.l) 
* a2(X)/ 

The 'first level algebraic Bethe ansatz' will be the same in both cases. The difference 
comes in only on the second level. 

The first step of our calculation is to rewrite the graded Yang-Baxter algebra ® 
with i?-matrix |_-) (X) in block form: We introduce the shorthand notations 

5(X) = (5i(^),52W), C{1) = 
Then the 3x3 monodromy matrix T (k) can be written as 

The defining relations of the graded Yang-Baxter algebra ® can be read as a 9 x 9 
matrix equation. Let us denote the n x n unit matrix by /„. A similarity transformation 
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with the matrix 



X 



(h 



\ 




\ 



hi 



(B.4) 



which cyclically permutes the 5th, 6th and 7th row and column, followed by a simi- 
larity transformations which multiplies the 5th, 8th and 9th row and column by — 1 , 
transforms this 9x9 equation into 



/I 



bh ah 
ah bh 



\ 



(Ai 
Ai 
C( 

\Ci 



Ac 
C( 



A®giB 



C AC 
fA Cc 



5( 



B D( 



\C®g>C C®giD D( 



'A 
C 
'A 
iC 



^lA Ac 

^'C Ac 

g/A Cc 

t,i C C C 



B 



B®gfA 



D 5C 
fB D( 



C 



B( 
5C 
A D( 
C DC 



5C 
5C 
D( 
DC 



)g'B\ 
)g,D 



/I 



V 



bh ah 
ah bh 



\ 



(B.5) 



For the formula to fit on the line we suppressed the arguments and adopted the follow- 
ing convention: X = X{X), X = X{ij) for X = A, . . . ,D. Moreover, a = a{iJ — X) and 
b = b{n-X). The 4x4 matrix 



/I 



v 



b 
a 



a 
b 



\ 



b — aj 



(B.6) 



satisfies the Yang-Baxter equation. It is related to a special case (rank 2, grading 
— )) of the i?-matrix ([T]) by equation ^ and is therefore unitary. 



r(X)r(-;i) =/4. 



(B.7) 



The grading {+, — ), corresponding to p{\) =0, p{2) = 1, also appears in the graded 
tensor products in ( IB.5I ). where it is abbreviated as g'. These graded tensor products 
are defined by equation (IA.4I) which makes not only sense for square matrices, but for 
arbitrary nxm matrices. Thus, thinking of B{X) as a 1 x 2 matrix with row index 1, and 
similarly of C{X) as a 2 x 1 matrix with column index 1 and of A (X) as a 1 x 1 matrix 
with row and column index 1, all graded tensor products in (IB. 51 ) are well defined. We 
have, for instance. 



B,il)C\^) B2{X)C\^i) 
B,{X)C\^) -B2{l)C\ij) 



(B.8) 
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We would like to remark that the defining relations of the graded Yang-Baxter algebra 
of the gl(l|2) model, when written in block form (IB.5I) . resemble the corresponding 
relations for the gl(2) model. 

Out of the 16 relations contained in (IB.5I) we shall need the following 4 for the first 
level algebraic Bethe ansatz, 

B{X) ®g^B{p) = {B{ij) (^g'B{X))r{p-X), 



(B.9) 
(B.IO) 



(B.ll) 
(B.12) 



a{n — X) a{i-i — X) 

f{p-l){D{l) ^g'D{p)) = (D(^) ®g,D(X))r(A/-^). 

Note that, by (IB.9I) . B{X) constitutes a representation of the Zamolodchikov algebra, 
and, by (IB. 121) . D{X) is a representation of the Yang-Baxter algebra of the gl(l|l) 
model. 

Our goal is to calculate the eigenvectors of the transfer matrix t{X) = —A{X) — 
strg/(D(?i)). In analogy with the gl(2) case we shall first of all consider the com- 
mutation relations of a multiple tensor product B{Xi) ®g' ■ ■ • (S>g' B{Xm) with A{X) and 
sir gi {D(k) ) . These commutation relations can be obtained by iterating equations (IB. 101) 
and dBlTT) : 

N N N 

(S)g>B{lj)]A{l)Yl 

■7=1 

A' 



A{1) lS),B{Xj) 

o 

N r 

N 



1 



fJ-^aiX-Xj 
^,B{Xk) \Sj-iAiXj 



k=l 



biX-Xj' ^ 



n 



1 



a{X-Xj) f^J[ a{Xj-Xk) 

k^j 



(B.13) 



r ^ if r ^ 1 1 ^ 1 

7=1 7=1 ■'^^ ^ 



N 



N 



■Y,{h®g'B{X)®g: ® ,B{Xk) \PQi{h®g'Sj-i)- 

k^j 

b{Xj-X) ^ 



{h®g:^iXg:{T{Xj))] 



n 



1 



a{Xj-X) j^J^ a{Xk-Xj 

k^j 



(B.14) 



Here the operators B{X) in the multiple tensor products are multiplied in ascending 
order. T {X) is defined in equation (l29b . The operators 5j_i appearing on the right 
hand side of (IB. 131) are given as 

5,._i = [fiXj - Xi) 0g, /f ^-')) . . . (/f ^'-2) rih - ^7-1) ^s' if'"''^) , (B.15) 
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for j = 2, . . . ,N. We further define So = id. The use of the graded tensor product 
in (IB.15I) makes sense, since all non-zero matrix elements appear in such a way that 
they can be interpreted as even elements of End(i^). Pqi is a graded transposition 

operator [36] on (End(C^))®''^^^'' regarded as a graded algebra (see appendix IXI). In 

the canonical basis {eo^' , . . . , e^a^^} it is expressed as 



(B.16) 



with 7 = and k = I. In general the operators Pjk induce the action of the symmetric 

group on the site indices of the basis elements e/^. 

Equations (IB.13I) and (IB.14I) can be proven by induction over (see [1]). As 
compared to the case g = (+,—,—) treated in [1] the main difference and the main 
difficulty was to properly treat the graded tensor products with the induced grading 
g' in the derivation of equation (IB.14I) . The remaining part of the derivation is now 
similar to the corresponding calculations for the case g = (+,—,—). We shall discuss 
it only briefly. 

We take the super trace in space zero of equation (IB. 141) and subtract it from minus 
one times equation (IB. 131) . Taking into account that = — A(?i) — strgf{D{X)) for 
the grading g = (—,—,+) under consideration we obtain 



N 



N 



A' 



7=1 



I) 



7=1 



N 



k=l 



1 



k^i 



N 



k+j 



(B.17) 



This is a form rather typical of an algebraic Bethe ansatz calculation with wanted and 
unwanted terms on the right hand side of the equation. 

The operator str^/ (T (X) ) acts on the tensor product (C^) '^^ ®J-[ . Its eigenvectors 
are independent of X since T iX) is a representation of the Yang-Baxter algebra with 
7?-matrix f{X), 



(B.18) 



This holds first of all for TgtQO) by construction (see (OTl) ) and for D{X) by equation 
(IB. 121) . It still holds after inserting the factors (o^)®^ into D{X) which then becomes 
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© (X) (see equation (|30l)). But, due to the factors (a^)®^ the entries of -D (k) and Tgi (k) 
super-commute, and (IB.18I) holds because of (IA.5I) . 

Following Kulish and Reshetikhin [11] we define the 'vacuum subspace' J-[q d 9{ 
by the conditions 

A(X)4) = ai(X)4>, (B.19a) 
C(X)4> = 0, (B.19b) 

for all 4> G :?/o- Clearly, J-[q is a linear subspace of H . The following lemma [12] can 
be proven in a similar manner as in [1]. 

Lemma 1. O-Cq is invariant under the action ofD(k). 

Equivalently we may say that the space spanned by all linear combinations of vec- 
tors of the form D\{fii) . . .D\{iim) is a linear subspace of ^Hq. 

Suppose (pt2 G (C^)®^(S)i^o is an eigenvector''! of str^/('T' with eigenvalue 
A(?i). Then cpQ is a 2^-column vector with entries in 9{q, and (IB.19al) holds for this 
vector, 

A(;i)cpa = ai(X)(pf2. (B.20) 

Thus, cpi^ is an eigenvector of the operators in curly brackets on the right hand side 
of (IB.17I) . Since the graded tensor products of vectors B{'kj) form a 2^-row vector we 
conclude that 

A' 2 

(g)^,5(?i;)](pt2= £ B,^{U)...Bi^{kMW-'''^&^ (B.21) 
is an eigenvector of the transfer matrix if the Bethe ansatz equations 

^=1 a[kj-kk) 

are satisfied, which is just the condition for the unwanted terms in the second curly 
bracket on the right hand side of ( IB. 171 ) to vanish. The corresponding eigenvalue of 
t{k) is 

A{k) = -a,{k) n - A(^) n -71^, ■ (B-23) 

fj^a{k-kj) fj^a{kj-k) 

The remaining task is to solve the eigenvalue problem of strg/('T (k)) on the space 
{C^)^^^o. This task can be accomplished by a second Bethe ansatz which is 
possible, because T {k) is a representation of the Yang-Baxter algebra (see (IB.18I) ) and 
the vector 

a=(y^f2 (B.24) 

""We loosely write cpil instead of (p i2. 
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is a highest vector for T (k) . In fact, introducing the explicit form of the L-matrix (l33t 

and the explicit form of |_), equation (I30cl) . into the definition (l29l) of "T (X) we 

obtain 



(B.25) 



For the construction of the eigenvectors of str^/('T (X)) we extract the following 
commutation relations from (IB.18I) . 



(B.26a) 



(m - mm = _ {am - m) ■ (B.26b) 



a{X — i-i) 



Here we referred back to the notation for the matrix elements introduced in (l29l) . Iter- 
ating ( IB.26bl) we obtain 



M 

str.K'TlX)) ]\B{^ik) 
M , 



M 



. M -, Ml 



k=l 



1=1 
l^k 



1 



(B.27) 



where the products over the B{^k) are ordered in ascending order and by definition 



k-\ 



1=1 



for k = 2, . . . ,M and 5o = 1 • It follows that 



(B.28) 



(B.29) 



is an eigenvector of str^/(T (k)) with eigenvalue 



A' -.Ml 

A{1)= a2{l)-a3{l)Yl<^J-^) riTnT^ 



if the Bethe ansatz equations 

asiPk) 



N 

n 



k=l,...,M, 



(B.30) 



(B.31) 



are satisfied. Inserting the expressions ( IB. 291 ) into (IB.21I) and (IB. 301 ) into ( IB.22I) . (IB.23I ) 
we arrive at the results shown in sectionHJand our derivation for the case g = (—,—, +) 
is complete. 
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It is now relatively easy to perform the algebraic Bethe ansatz for the remaining 
grading (—,+,—). According to our remark at the beginning of this appendix this 
case is equivalent to considering the grading (—,—,+) with the vacuum action shown 
in equation dB.ll) . Therefore only the second level Bethe ansatz, starting from equation 
(IB .241) has to be modified. Since D{X) now acts as a lower triangular matrix on Q. we 

must choose the auxiliary vacuum for the second Bethe ansatz as t2 = ( J Q. Then 
T (X) has lower triangular action on Cl, and an algebraic Bethe ansatz (with C replacing 
B) becomes again possible: Introducing the explicit form of the L-matrix (l33t and the 
explicit form of ^ ), equation (I30bl) . into the definition of T (X) we obtain 





yv ajlj-X) 



a. 



(B.32) 



= 1 a{X~Xj) 



where we used the identity a(^) —bCk) = a{'k)/a{—X). Instead of the commutation 
relations (IB. 261) we now need the commutation relation between C{X) and C(/j) and 
between strg/{T {X}) and C(/j). These commutation relations are again contained in 
(IB. 181 ). They are of the same form as (IB. 261) with C replacing B and the arguments X, 
}i of the functions a and b interchanged. It follows that 



(B.33) 



is an eigenvector of str^/(T (k)) with eigenvalue 



A{X)= a3{X)Yla{Xj-X)-a2{X)Yl- .WW 

L fj^ fj-^ aiX - Xj) J J^J^ a{pk - X) 



n 



1 



if the Bethe ansatz equations 



N 

n 



k=l,...,M, 



(B.34) 



(B.35) 



are satisfied. Inserting the expressions (IB.33I) into (IB. 211) and (IB. 341) into (IB.22I) . (IB. 231) 
we obtain the results for ^ = (—,+,—) shown in section IH 
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